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ABSTRACT 



A new approach to defining inlet distortion factors for axial -flow 
compressors is presented. A summary of past analyses of inlet distortion, 
along with a brief review of flow-distortion definition and compressor 
performance, provides a background for the proposed theory. Inlet flow- 
distortion in terms of total -pressure fluctuations at the compressor face 
is converted to vorticity. The effects of vorticity on the internal 
aerodynamics of the compressor are investigated with some approximate 
calculations included to indicate the validity of this approach. Refine- 
ments and alternate solutions to the theory which could lead to definition 
of a more comprehensive and reliable definition of stall-inducing inlet 
distortion are included. 



2 



TABLE OF CONTENTS 



I. INTRODUCTION — 10 

II. STATEMENT OF PROBLEM 12 

A. AXIAL COMPRESSOR PERFORMANCE 12 

B. PRESENT STATE OF INLET DISTORTION ANALYSES 18 

C. VORTICITY APPROACH 21 

1. Compressor Stall 21 

2. Secondary Flows in Axial Compressors 23 

3. The Theory of Vorticity Analysis 25 

HI. CONVERSION OF PRESSURE TO VORTICITY MAPS 28 

A. CROCCO'S THEOREM — 28 

B. ORDER OF MAGNITUDE ANALYSIS 31 

C. APPLICATION OF RESULTS 39 

IV. ANALYSIS OF VORTICITY MAPS 46 

A. RADIAL VORTICITY INLET EFFECTS 46 

B. CIRCUMFERENTIAL VORTICITY INLET EFFECTS 51 

C. PASSAGE AND EXIT EFFECTS 56 

D. CALCULATION OF SECONDARY CIRCULATIONS 61 

V. SUMMARY AND REVIEW - 72 

VI. CONCLUSIONS 75 

APPENDIX A ORDER OF MAGNITUDE ANALYSIS 77 

APPENDIX B CONVERSION OF Pj MAP TO VORTICITY - J85 96 

APPENDIX C ESTIMATION OF SOME J85 PARAMETERS Ill 

APPENDIX D CALCULATION OF ABSOLUTE VORTICITY 116 

APPENDIX E CALCULATION OF CIRCULATION FROM Wq MAP FOR J85 120 



3 



LIST OF REFERENCES - — 125 

INITIAL DISTRIBUTION LIST 131 

FORM DD 1473 - 138 



\ 



4 



LIST OF TABLES 



I Characteristic Times 30 

II Relative Orders of Magnitude 37 

III Possible Combinations for Radial Inlet Vorticity 46 

IV Possible Combinations for Circumferential Inlet Vorticity 52 

V Magnitudes of Circulation Components 71 

B-I , Pressure Data Taken in Radial Direction from 

J85 Total-Pressure Map 98 

B-II Pressure Gradients Obtained from Slopes in Fig. B-2 102 

B-III Pressure Data Taken in Circumferential Direction 

from J85 Total-Pressure Map 104 

B-IV Pressure Gradients Obtained from Slopes in Fig. B-3 107 

B-V Net Fluxes of Vorticity in Three Coordinate Directions 110 

D-I Absolute Vorticity for Three Circumferential 

Positions on ug Map 118 

E-I Relative Vorticity and Circulation Components at Radial I 122 

E-II Relative Vorticity and Circulation Components at Radial II — 123 

E-III Relative Vorticity and Circulation Components at Radial III — 124 






5 










$>• 




LIST OF FIGURES 



1. Compressor Performance Map, J85-6E-13 

2. Compressor/Receiver Matching Characteristics 

3. Radial Pressure Distortion 

4. Circumferential Distortion 

5. Loss in Steady-State Distortion Tolerance with Turbulence — 

6. Typical Stall Patterns 

7. Stall Cell Growth and Rotation in a Blade Row 

8. Secondary Flow and Vortices in an Axial -Flow 

Compressor Rotor 

9. Nomenclature 

10. Radial Vorticity Contours, J85 

11. Circumferential Vorticity Contours, 085 

12. Contours of Net Flux of cd'^ for 085 

13. Radial Vorticity Contours, TF-30 

14. Circumferential Vorticity Contours, TF-30 

15. Radial Vortex-Sheet Passing Through a Stator 

16. Radial Vortex-Sheet Passing Through a Rotor 

17. Radial Vortex Filaments Passing Through a Stator 

18. Radial Vortex Filaments Passing Through a Rotor 

19. Circumferential Vortex-Sheet Approaching a Stator 

20. Components of Circumferential Vortex-Sheet 

Approaching a Rotor 

21. Streamwise Component of a Circumferential Vortex-Sheet 

at a Rotor 

22. Circumferential Vortex Filament Approaching a Stator 

23. Vortex Filament Passage Through a Cascade 



12 

13 

15 

16 

19 

22 

22 

25 

28 

41 

42 

43 

44 

45 

47 

48 

50 

51 

52 

53 

54 

55 

57 



6 



24. Secondary Flow Downstream of a Blade Row 61 

^5. Nomenclature 62 

26. Vector Transformation From Absolute to Relative Vorticity 64 

27. Passage of a Vortex Filament Through a Rotating Cascade 65 

28. Data Points for Calculation of 66 

29. Vorticity Components at Radial I. ^ 67 

30. Exit Secondary Circulation Components at Radial I 69 

31. Flow Diagram for Analysis; Steady Flow 73 

32. Composite Theory for Analyzing Vorticity Maps 74 

A-1 Peak Instantaneous Pressure Contours, J85 78 

A-2 Transformed Pressure Contours 80 

A-3 Frame One of Instantaneous Pressure Contours Sequence, TF-30 -- 81 

A-4 Frame Two of Instantaneous Pressure Contours Sequence, TF-30 — 82 

A-5 Frame Three of Instantaneous Pressure 

Contours Sequence, TF-30 83 

A-6 Frame Four of Instantaneous Pressure 

Contours Sequence, TF-30 84 

A-7 Frame Five of Instantaneous Pressure 

Contours Sequence, TF-30 85 

A-8 Frame Six of Instantaneous Pressure 

Contours Sequence, TF-30 86 

A-9 Frame Seven of Instantaneous Pressure 

Contours Sequence, TF-30 87 

A-10 Frame Eight of Instantaneous Pressure 

Contours Sequence, TF-30 88 

A-11 Frame Nine of Instantaneous Pressure 

Contours Sequence, TF-30 89 

A-1 2 Sources of Dissipation Within Inlets 90 

A-1 3 Pressure History for Steady Inlet Flow 91 

A-1 4 Shock Wave Twitter 91 

A-15 Total-Pressure Fluctuations From Oscillating Shock Waves 92 



7 



A-16 Pressure Waves in the x-t Plane as a Result of 



Shock Wave Twitter 93 

A-17 Nomenclature 94 

B-1 Pressure Data and Scales for J85 97 

B-2 P' vs. R' 99 

B-3 P' vs. 0 105 

B-4 Method of Divergence Application 108 

C-1 Blade Passage Dimensions 112 

C-2 Velocities Ahead of Rotor 114 



8 



ACKNOWLEDGEMENTS 



The author sincerely appreciates the patient and unselfish guidance 
offered by Professor Allen E. Fuhs of the Department of Aeronautics at 
the Naval Postgraduate School during the completion of this thesis. 

The author also gratefully acknowledges the diligent work performed by 
his wife, Sally. 

Several persons and organizations supplied data for this work. The 
author wishes to thank Mr. Robert P. Zalis, Mr. Paul H. Kutschenreuter, 
Jr., and Mr. David Jamison all of General Electric Company, Mr. Forrest 
Schubert at the AFAPL/WPAFB, and Mr. James Lonsdale of Pratt and Whitney 
Aircraft for their willingness in providing data and information. 

Without their help this thesis would not have been possible. The use 
and interpretation of the data are, of course, the responsibility of the 
author. 



9 



I. INTRODUCTION 



Engine stall has always been a problem of varying magnitude in jet 
powered aircraft employing axial -flow compressors. At present aircraft 
gas turbine performance is largely limited by flow separation and surge 
in the compressor. Generally, highest compressor efficiencies occur 
near the point of flow breakdown (stall), and so the design margin 
betvjeen normal operation and stall is kept to a practical minimum. A 
compressor's ability to remain below this stall limit is significantly 
affected by the amount and degree of distortion in its inlet flow. 

Inlet flow distortion has been called by many names, all of which basi- 
cally describe variations in total -pressure , total -temperature, or 
velocity over the inlet cross-section. Previously, inlet flow distor- 
tion has been at a level which the compressor could handle without a 
serious degradation in performance. However, with the greater speeds 
and performance demanded of today's high-specific-energy turbomachines, 
inlet distortion at the compressor face can alone be great enough to 
cause compressor stall and surge. One prominent source of distortion 
has been traced to shock-wave/boundary-layer interactions within the 

inlet. Others are steam ingestion, pressure transients due to afterburner 

\ 

light-off, ordnance discharges, etc. The magnitude of the problem in 
some recent aircraft designs has caused much research to be conducted in 
this area, specifically correlating maps of total-pressure variations at 
the compressor face to engine surge. The results of this work have 
yielded an immense amount of empirical data. However, reduction of these 
data has yet to provide a technique which adequately describes the distor- 
tion accommodation of a given engine-inlet configuration over its entire 
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range of operation. In contrast to empirical correlations yielding a 
statistical yes-no answer to the question of whether or not the engine 
will stall, methods are needed which link stall criteria for a blade 
element to distortion patterns. This paper presents a new approach to 
analyzing already existing data and provides evidence of a better 
insight into the actual fluid mechanics involved in the problem. 
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II. STATEMENT OF PROBLEM 



A. AXIAL COMPRESSOR PERFORMANCE^ 



The performance of an axial-flow compressor such as found on air- 
craft turbojet engines is best described by its steady-state performance 
map. Figure 1 illustrates a typical compressor map. "Pressure ratio" 
is exit total -pressure to inlet total-pressure and "corrected airflow" 
is engine mass-flow rate corrected for total-temperature and total- 



pressure. 



O 

I- 

< 

a: 




FIG. I COMPRESSOR PERFORMANCE MAP 
J85- GE- 13 [adapted from Ref.z] 



^The information and examples in this section were adapted from a 
report by Brimelow [Ref. 1]. 
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The surge line defines the maximum pressure ratios which the compressor 
can achieve throughout its range of rotor speeds. Operation above this 
line is accompanied by a drastic reduction in pressure ratio. This is 
called "surge," and it is caused by excessive aerodynamic blade stall. 

The stability of a system is determined by the slope of its pressure 
ratio/flow characteristic. Consider steady operation at point M-j of 
Fig. 2. 




FIG. 2 COMPRESSOR / RECEIVER 

MATCHING CHARACTERISTICS 

Now if the flow became unsteady and a perturbation were to cause a 
reduction in airflow to point S-j (constant speed), the compressor would 
be delivering an instantaneous pressure greater than that in the 
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receiver. Effectively the compressor back-pressure has been lowered, 
and the flow would increase until a match is again reached. Similarly 
for an increase in flow to point $2, the increased back-pressure would 
cause a return to stable operation at point M^. However, if the compres- 
sor is matched on the positive slope of the characteristic, by the same 
reasoning, unsteady conditions exist and perturbations are aggravated 
to "surge." 

A stall margin, defined as the vertical distance between the nominal 
(design) operating line and the surge line, can be thought of as a 
measure of the compressor's ability to accommodate flow transients which 
instantaneously load the machine. This margin is dictated by overall 
engine-inlet design and is affected by many different phenomena. Unfor- 
tunately the point of maximum efficiency in a compressor is just prior 
to stall or just below the surge line on the performance map, and there- 
fore a trade-off between efficiency and surge margin is required. The 
expected flow transients which the compressor must be able to accommodate 
will affect the design efficiency. 

Flow transients initiated within the engine system are normally 
caused by power level (throttle) movement. Transients can also be 
initiated by flow distortion at the compressor face. Of the various 
flow distortion parameters, Brimelow [Ref. l] states that inlet total- 
pressure and total -temperature have been identified as having the greatest 
effect on engine stability, which is to a large extent the compressor's 
ability to remain below the surge line. Total-pressure variations can 
be steady-state spatial distortion, time-varying spatial distortion, 
in-phase pulsations which are spatially uniform over the compressor face, 
or a combination of these. Low frequency spatially uniform (planar) 
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pressure fluctuations (on the order of 20 CPS) do not affect the com- 
pressor but do have a detrimental effect on the inlet and engine control 
systems. High frequency planar fluctuations (more than one per revolu- 
tion) affect compressor stability much like time-varying spatial 
distortion. 

Spatial distortions can be further segregated into radial and cir- 
cumferential inlet flow distortion. This separation into circumferential 
and radial distortion was arbitrary, somewhat motivated by geometric 
simplicity. A radial total -pressure distortion such as shown in Fig. 

3(a) produces very little variation in inlet static pressure. The net 
result is a low velocity in the low pressure region and a high velocity 
in the high pressure region. 




TIP 



LOW PRESSURE 



HUB 



TIP 



(a) 



(b) 



FIG. 3 



RADIAL PRESSURE DISTORTION 
from Ref. l] 
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In the low velocity region the blades will operate at higher incidences 
and higher total-pressure ratios. The net result will be a reduction in 
the compressor stall margin. 

For the circumferential type distortion shown in Fig. 4(a), the 
associated surge line reduction is a function of the basic stage-match 
of the compressor, the axial spacing between blade and vane rows, and 
the unsteady response of the rotor blades. Since the blade rows restrict 
circumferential redistribution of flow, each circumferential segment 
performs essentially independently. 






q 



Undistorted 

Inlet Flow 



Low Pressure 



(a) 



180® 

(b) 



360® 



FIG. 4 



CIRCUMFERENTIAL DISTORTION 
from Ref. l] 




All segments are required to deliver the same exit static pressure so 
that in the area of low inlet pressure there is a high pressure ratio 
and vice-versa. Correspondingly, the high pressure ratio forces a low 



axial velocity in this area and a high velocity in the segment of high 
inlet pressure. These variations in velocity, produced less than one 
duct-diameter upstream of the compressor, result in static pressure 
fluctuations at the compressor face. This mechanism of varying veloci- 
ties (incidences) can induce surge when the average pressure ratio of 
the distorted flow is less than the surge pressure ratio for undistorted 
inlet flow. The effects of circumferential pressure pulses on a TF-30 
engine are presented in Ref. 3. It was found here that the amplitude of 
a pulse necessary to stall the engine varied inversely with pulse dura- 
tion and that the engine was most sensitive to distorted sectors of 180 
to 240 degrees. 

The distortions which are encountered in actual compressor applica- 
tions such as jet aircraft inlets are neither radial nor circumferential 
but a combination of the two. Furthermore the magnitude, location, 
shape, and duration of these pressure variations change continually. 

It is virtually impossible to evaluate precisely the effect that these 
distortions have on engine stability because there are so many distor- 
tion patterns significant to each compressor surge. Therefore it has 
been necessary through experimental programs to develop generalized 
correlations between pattern characteristics and corresponding effects 
on compressor stability. Correlating the shape, location, spatial 
extent, and magnitude of the inlet distortion requires analyzing many 
radial and many circumferential patterns and relating the overall picture 
to engine stability. From this is generated a distortion factor which is 
a measure of the associated loss in stability-margin for that given 
propulsion system. Because compressor geometry and blading are not 
considered in arriving at this parameter, the result cannot be applied 

I 

to other systems. 
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For a more detailed and complete analysis of this subject, consult 
Brimelow's excellent work in Ref. 1. 

B. PRESENT STATE OF INLET DISTORTION ANALYSES 

Until a few years ago, engine-inlet compatibility limits were success- 
fully determined by evaluating an engine's tolerance to steady-state 
distortion. Normally this was done with screens or similar obstructions 
placed in the inlet and the compressor face instrumented with low-response 
pressure probes. The engine was operated under all conditions of interest, 
and by comparing pressure data at the compressor face with the engine's 
performance, distortion factors were derived. It was known that unsteady 
flow was being produced by these test devices, but no attempt was made to 
evaluate its effect. However, when this technique was applied to the 
TF-30 turbofan engine during the F-lllA flight test program, the data did 
not correlate well nor was the method found valid for other fan type 
engines at that time. An improved steady-state distortion factor, 
was developed, but it too fell short of providing the needed total evalu- 
ation of inlet distortion. It was becoming more evident that there was 
another significant variable involved, namely turbulence. Turbulence can 
be described as flow with random pressure fluctuations and can be quanti- 
fied as a fluctuation in RMS total -pressure , (aP|^ 2 )/(Pj) • With 
increased aircraft speeds, turbulence has taken on greater importance 
since there is a significant amount of turbulent flow associated with 
the boundary-layer/shock-wave interaction within supersonic inlets. It 
was felt that steady-state distortion indexes were still of value, but 
a correlation with turbulence was necessary. One such empirical relation 
showing a decrease in distortion tolerance, increasing turbulent 

energy is depicted in Fig. 5. 
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TURBULENCE - (AP^^s 

FIG. 5 LOSS IN STEADY- STATE DISTORTION 
TOLERANCE WITH TURBULENCE 
[adapted from Ref. 4 ] 

Turbulence has two primary effects on axial compressor performance. 
First, turbulent inlet flow causes a reduction in compressor surge line, 
its sensitivity being a function of turbulent frequency and amplitude. 
Second, a reduction in flow capacity accompanies inlet turbulence. These 
two phenomena were observed in tests performed by Brimelow and Plourde 
[Ref. 5] on a TF-30 turbofan engine. This all led to the development of 
new distortion indexes such as the and parameters of Ref. 6. 

However, lack of correlation caused Burcham and Hughes to raise doubt as 
to the usefulness of the conventional turbulence definition, • 

Bellman and Hughes [Ref. 7] found no consistent relationship at all 
between turbulence factors and surge, and in some cases surge was experi- 
enced as turbulence level was decreasing. Therefore it was concluded here 
that peak-to-peak pressure fluctuations at key areas of the compressor 
face are more important in surge prediction than are turbulence levels. 
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Also extensive experimental data have been obtained from a J85-6E-13 
turbojet placed in a supersonic wind-tunnel running at Mach numbers 
close to 2.5. By varying the characteristics of the engine inlet duct, 
a desired level of time variant flow distortion was produced through 
internal shock interactions [Ref. 2]. Calogeras, Burstadt and Coltrin 
found for this J85 configuration that time-averaged distortion could be 
well below a stall-inducing value while instantaneous peaks exceeded 
this value. In addition, severe pressure fluctuations could result in 
an instantaneous distortion of a type quite different from the time- 
averaged distortion. 

Of these few representative programs discussed, there are two 
striking similarities worthy of note: (1) Each is dedicated to using 

inlet flow distortion described by raw total -pressure data to predict 
compressor surge. No attempt at determining the cause or nature of the 
stall is made. (2) There is much uncertainty and disagreement associated 
with the present methods of predicting distortion induced stall. A 
universal distortion parameter has not been agreed upon. 

Obviously a better understanding of the internal fluid dynamics, 
while a difficult task at best, would nonetheless allow a more compre- 
hensive generation of a stall-inducing distortion parameter. Total- 
pressure is the easiest and perhaps the only flow parameter which can be 
measured at the compressor face with sufficient accuracy and frequency 
response to gain meaningful data. The vast amount of pressure data 
available due to the previously mentioned programs provides an exten- 
sive empirical data base with which to treat this problem. In 
evaluating these data, inclusion of time and temperature dependence 
could only improve the usefulness of the results. Therefore it is 
proposed that vorticity become the basis of a "universal distortion 
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factor." This thesis will show how existing total-pressure maps along 
with total -temperature and time data can be transformed into vorticity 
maps. Analysis of the resulting vorticity maps will provide evidence 
as to the flexibility of this method. 

Vorticity is more closely related to fluid dynamics than is total 
pressure alone. Compressor stall and engine surge are flow-separation 
phenomena; therefore the key to understanding their causes lies in the 
internal dynamical situation. 

C. VORTICITY APPROACH 

1 . Compressor Stall 

The most significant phenomenon affecting the performance of an 
axial-flow compressor is stall. If the flow rate through a compressor 
blade row is decreased while the rotor speed is maintained constant, 
angles of attack on the blades will increase and eventually flow sepa- 
ration will occur. This condition is called stall and is similar in 
some respects to the stall of an isolated airfoil; differences are due 
to the adjacent blading in a cascade. Stalled operation of a compressor 
is accompanied by a drastic decrease in efficiency and excessive oscil- 
lating blade loads capable of causing structural failure. Stall occurs 
in the range of unstable compressor operation and hence can trigger 
surge. Usually a compressor stalls in an asymmetric manner; that is, 
well defined regions of the annulus are stalled while the remaining area 
is unstalled; see Fig. 6. The flow through these stalled patches is 
severely restricted in accordance with the decrease in total flow-rate 
through the compressor, thus the flow is diverted around this blockage. 
Effectively this increases blade angle-of-attack on one side of the 
restriction while decreasing it on the other side. In Fig. 7 blade No. 1 
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FIG. 6 TYPICAL STALL PATTERNS 









Decreasing 
Angle- of - Attack 



Increasing 
Angle- of - Attack 



FIG. 7 STALL CELL GROWTH AND 
ROTATION IN A BLADE ROW 
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•will soon become unstalled while at blade No. 5 flow separation will 
commence. This condition is termed "rotating stall" since it propagates 
circumferentially in a direction opposite to rotor rotation. In an 
absolute frame of reference this speed of circumferential propagation 
is 30-60 per cent of the rotor speed [Ref. 8]. When these conditions 
become severe enough to cause net flow-rate fluctuations with time, the 
situation is called "surge." This describes operation in the unsteady 
area of the compressor performance map in Fig. 1. Valensi [Ref. 9] has 
made a thorough investigation into these phenomena with some excellent 
photographic results. Of the various analytical methods available for 
predicting compressor stall, the most generally accepted employs the 
limiting-blade-loading criteria of Lieblein, Schwenk, and Broderick 
[Ref. 10]. Their "limiting diffusion factor," D, is a function of the 
velocity relative to the blade and the circulation about the element. 

If a constant axial velocity assumption is made, this simplifies to a 
function involving the flow turning angle and cascade solidity. For 
design purposes where such stage parameters may not be available, 

Lieblein [Ref. 11] has derived an equivalent diffusion ratio, which 
is based on inlet and outlet conditions to the compressor and the type 
of blading installed. Unlike a single airfoil where angle-of-attack is 
used to predict stall, the critical variable for a cascade of blades is 
the flow turning-angle. The exit-flow angle as well as the inlet angle 
must now be considered. 

2. Secondary Flows in Axial Compressors 

Returning now to normal compressor operation, flow patterns 
based on axisymmetric flow distribution and an infinite number of blades 
have been described quite well in the literature. However, it is necessary 



23 



to consider a finite number of blades since there are wall boundary- 
layer interactions, tip clearance leakages, blade boundary-layers, and 
centrifugal effects which are not represente°d adequately by assuming 
axisymmetric flow. Although these asymmetric phenomena have been 
studied independently for reasons of simplicity, it must be remembered 
that their effects are not necessarily additive (non-linear) and may 
have strong interactions. 

Figure 8 illustrates the various secondary flows which have 
been found to exist in axial-flow turbomachines as a result of the above 
mentioned viscous and three-dimensional effects. Lakshminarayana and 
Horlock [Ref. 12] have made athorough review of secondary flows in axial- 
flow turbomachines and have indicated where in the literature more 
specific information can be found. Woods [Ref. 13] presents the current 
state of the art in analytical treatment of secondary flows and associ- 
ated losses in axial-flow turbomachines. 

The induced velocities caused by secondary flows are considered 
losses since only a small part of their kinetic energy is recoverable 
downstream. Secondary velocities may perturb the exit flow angle from 
a blade row. Near annulus walls where lies the origin of the central 
vortex in Fig. 8, the induced exit flow angle can cause local blade 
stall. Much work has been performed in analyzing these phenomena, and 
some analytical methods for predicting three-dimensional exit flow 
angles have been proved valid through experimental results [Refs. 14 
and 15]. 

Another source of strong secondary flows in blade passages lies 
in non-uniform inlet flow such as described previously. Treating inlet 
flow distortion as vorticity provides a natural basis from which to 
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FIG. 8 SECONDARY FLOW AND VORTICES IN AN AXIAL 
FLOW COMPRESSOR ROTOR [adaptedfromRef.il] 

compute resultant secondary circulation at the exit plane. The effects 
of this secondary flow on cascade performance can be studied analogously 
to boundary-layer induced circulation. Comparison with an appropriate 
blade-loading criterion may predict the compressor's unstable operation. 

3. The Theory of Vorticity Analysis 

The thermodynamic properties v/hich are significant in compressor 
performance are total -pressure (velocity), total -temperature (enthalpy), 
specific heat, and molecular weight. Variation of any one or combination 
of these properties with time and/or position at the compressor face 
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constitutes inlet flow distortion and should be an integral part of any 
universal distortion parameter. To date distortion factors dealing only 
with one inlet property at a time have been reported, i.e. total - 
temperature variation ->■ Kj ^ K-j-(stall) or total -pressure fluctuation ^ 

Kp -»■ Kp(stall). This type of piecemeal approach, while providing satis- 
factory results for certain restricted applications, cannot be expected 
to remain valid under all conditions. In reality, an aircraft engine 
will encounter all of these fluid distortion properties at one time or 
another. A partial list of causes of these is: inlet characteristics, 

environmental conditions, rocket exhaust from ordnance, hot runway opera- 
tion, steam catapult launches, and the normal operating envelope of the 
aircraft. There is definitely need for a more realistic or universal 
distortion parameter, one which accounts for the total inlet distortion. 

It is the intent of this work to show that a distortion parameter based 
on vorticity would be more flexible and yield more comprehensive surge 
prediction criteria than present methods. 

Crocco's Theorem provides a versatile means of converting exist- 
ing inlet distortion maps to vorticity in that pressure, temperature, 
and time are all available as desired for inclusion in the transformation. 
As indicated previously, time is an important parameter in such an analy- 
sis, the aspects of steady-state versus turbulence induced stalls being 
quite distinct. Hawthorne [Refs. 16 and 17], Preston [Ref. 18], and 
Horlock [Ref. 15] all deal with the passage of vorticity through a cas- 
cade. Horlock and Hawthorne calculate induced exit flow angles which 
involves a solution to Poisson's Equation. Calculating appropriate exit 
circulation and comparing it with blade-lift circulation is a far simpler 
process that yields a qualitative feel for the relative size of this 

i 
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effect. Comparison of this secondary circulation with that caused by 
three-dimensional effects should also be made in order to look for 
possible interactions. 

The effect of inlet distortion on blade angle-of-attack must 
also be considered since this is the other "half" of the flow turning 
angle. This effect is much less subtle and can be dealt with in a more 
straight-forward manner. Areas of low pressure represent lower axial 
velocity, hence higher angles-of-attack. The unsteady response of the 
blades must be considered, particularly for circumferential inlet 
distortion where blade passing frequency is a factor. 



27 



III. CONVERSION OF PRESSURE TO VORTICITY MAPS 



A. CROCCO'S THEOREM 

In general form, Crocco's Theorem can be written as 



T grad S + u x curl u = grad H + -^ 

d t 



( 1 ) 



which holds in unsteady flow [Ref. 19]. From basic Thermodynamics 

Y - 1 \T^,/ 



S 2 - St / p 

— = - Ln 



R 



(-) 



TV 



Combine Eqns, (1) and (2), assuming H = Constant 



■RT ^ D 4. 1 

■p— grad Pj + u x curl u = -^ 
• T ' 9t 



(3) 



The following development will apply to the system defined in Fig. 9. 

r 




FIG. 9 NOMENCLATURE 
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Writing Eqn. (3) in component form using cylindrical coordinates yields 



-RT 8Pi 



Py 3r 



e + (u»oj - u_Un)e 
r ' e z z 9^ r 



8U 



r - 
- e. 



at 



(4) 



■RT 1 . 3U . 

Vz)4e - -eo 



(5) 



■RT 9Py^ yv 3U-,yv 

9 + (u,_h\. “ U ri\ )p = 9 

Pj ^ » t' Z — Z 



( 6 ) 



The following reference quantities taken at the compressor face will 
be used to non-dimensional ize the equations. Refer to Fig. 9. 



Py = P'Py 


Fy = Average Stagnation Pressure 


U = u'¥ 


a = Average Speed of Sound 


,D 

r = r 2 


D = Compressor Diameter 




0 ) = 0 )' 0 ) = Fluid Vorticity 

D 

Table I lists some of the characteristic times available. Blade passing 
frequency will be used in this development, since most available data 
include compressor speed; then 



t 



t' 



2ttR 




The numerical value of the reduced frequency for stall is based on the 
work of Carta [Ref. 20]. For values less than approximately 0.1, blade 
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TABLE I CHARACTERISTIC TIMES 



BLADE PASSING FREQUENCY 

V^ipN (1500) (36) 

2irR " ( 6 ) ( 2 ) 



= 4500 Hz 



ONCE-PER-REV 

V 

^TIP 

ZttR 



1500 
( 6 ) ( 2 ) "" 



WAVE PROPAGATION 



RADIAL — 

2R 


1200 

4 


300 Hz 




LONGITUDINAL 


a + u 


= 60,180 = 120 


Hz 


L 


STALL TIME 








REDUCED 


= 7T fC 


- 0.1 




FREQUENCY 


U 






f 


^(0.I)(I200) . 
(3) (0.2) 


Hz 
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stall is essentially equivalent to steady-state stall. For a reduced 
frequency larger than 0.1, stall is suppressed due to dynamic effects. 

Other work in this area has been performed by Carta [Ref. 22]. The 

$ 

paper by Ericsson and Reding [Ref. 21] is a survey of dynamic stall. 

Rewriting the component equations in non-dimensional form and 
canceling terms yields 



J_ 

yP' 



3p' 

3r' 







•) = 



1 3P' 

yPV 36 






1 3P' 

yP' 3z' 







(7) 

(8) 

(9) 



Equations (7), (8), and (9) are exact. Orders of magnitude of various 
terms will be examined to ascertain whether or not complexity can be 
reduced in a logical and consistent manner. 

B. ORDER OF MAGNITUDE ANALYSIS 
2 

Figure A-2 represents an inlet-flow distortion map at the face of 
a J85-GE-13 turbojet. This distortion map triggered surge. The map 
of Fig. A-3 is a similar situation in a TF-30 turbofan engine. These 



2 

Reference to figures and equations labeled with a letter/number 
combination have more detailed explanations which are contained in the 
Appendix. 
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maps will be considered representative of stall-inducing flow distortion 
for purposes of making an order-of-magnitude analysis on Eqns. (7), (8), 
and (9). 

Establish an order of magnitude, e , such that 




From Figs. A-2 and A-3 

- 0.2 1 e £ 0.2 

Between 12 and 5 o'clock. Fig. A-2, aP'/Sr' is a maximum. Refer to 
Appendix A for all numerical values. 

3P' 3(Pt/Pt) - .97 

= ! — — - = 1.53 

9r' 3(2r/D (7.1 - 6.05)(2/16.1 ) 

From Fig. A-3 at the 12:30 o'clock position 

3P' 1.04 - .96 

== = .685 

3r' (2.15 - 1.8)(2/6) 

3P' 

-0.7 < <1.5 

3r' 

Between 5 and 7 o'clock. Fig. A-2, a maximum value of 3P'/39 exists such 

that 1.20 - .95 

= = .416 per radian 

30 .6 rad 

Similarly from Fig. A-3 at 1 o'clock 
3P' 1.04 - .92 

= .343 per radian 

30 .35 rad 
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A discussion of the term ap'/az' is contained in Appendix A, Part II. 
An acoustic analysis taken at the compressor face yields Eqn. (A-2) 
which allows a magnitude estimation of aP'/9z'- 



aPf 



az 



face 




(A-2) 



ap' D ap.j. 

az' Ft az 




From Figs. A-3 through A-8 at the point "0", the following data are 
available. 



P' 

P' 



.92 

.96 



"T 

Pt 



P'Pt = 1480 psf 

1545 psf 



(Pj ■ P$) 



1/2 



p/2 



.0015/2 



= 475 ft/sec 



a = u/M = 475/. 4 = 1190 ft/sec 



aP' / 3 ft \ / 1 sec/ft \ /1545 - 1480 Ibs/fr \ 

az' \1610 lbs/ft2/ \475 + 1190/ \6 Fig x .2x10'^ sec/Fig/ 



ap' 

az' 



= .06 
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An average value of u^' for the J85 can be obtained from data in 
Appendix C. 

u 460 

u ' = ^ - = .415 

^ a 1110 



A maximum value would be 
P' = 1.17 = 



1/2 

(^T - ^S) 



p/2 



750 

u,' = = .675 

1110 



u 



I 

z 



= .7 



750 ft/sec 



Data with which to estimate the magnitudes of u^' and Ug' or their 
partial derivatives are not available. However, an upper bound of order 
e can be definitely established by considering the consequences of u^' 
or Uq' = O(e^). A radial velocity component of order unity is not 
possible within a duct. Any u^,' components present at the compressor 
face are caused solely by turbulent fluctuations and are probably of 
order however an upper bound of e will be assumed. A velocity 
component such as Ug', called pre-rotation, is often generated intention- 
ally as a means of establishing a desired angle-of-attack on the first 
blade row. The magnitude of this velocity should be smaller than u^' 
and certainly not larger. A velocity in the circumferential direction 
of order e^, about equal to the axial velocity, is unrealistic; and if 
it were present, the engine would just not run. For the case under 
consideration there is no pre-rotation so that any Ug' which does exist 

p 

is caused also by turbulent fluctuations and is probably of order e . 
Again an upper bound of e will be assumed. 
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Therefore 





Partial derivatives of u^' and u^' with respect to time can also be 
given only a bound with respect to u^'. 



There are two phenomena causing a time dependency. One is due to the 
motion of the rotor. For patterns of large extent on the Pj map, the 
appropriate frequency is once-per-rev. For patterns of limited extent, 
blade passing frequency is appropriate. Unsteady aspects of the flow 
cause a time dependency also. 

Continuing, the magnitude of cop' will be assessed. 



From Fig. A-2 between 8 and 9 o'clock, and data from Appendix C 



9u„' 9u ' 9u ' 



9t' 9t' 9t' 



(0 



r 




2a 



P' = .97 



u 




400 ft/sec 



z 



p/2 



P' = 1.17 



u 



z 



750 ft/sec 



A0 = 30° = .52 radians 



(350 ft/sec) (1.34 ft) 



(.52 rad)(.605 ft)2(1110 ft/sec) 



.814 per radian 



(0 



r 



.8 
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1 



Similarly, approximating Wq' yields 




From Fig. A-2 between 12 and 5 o'clock 



(JL). 




-(750 - 400) 
.06 ft 



ft/sec 



1.34 ft 

2(1110 ft/sec) 



-3.52 



0 ) 



0 



3,. 5 



For the axial component of vorticity 



(0 



Z 



2a 






D 

2a 



Since both partial derivatives are really unkown in size, this term 
will be evaluated later. 



The magnitude of au^'/^t' in Eqn. (9) can be estimated by again 
referring to the point "0" in Figs. A-3 through A-8. 




P' = .96 ~ ft/sec 

At = (6 figs.) X (.2x10"^ sec/fig.) ^ 

au^' TT /Au^x 3.14 / 374 - 231 

at' aRPSN\At/ (1190)(166) \ 1.2x10’^ 
au^' _ .06 
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Entropy can be expressed as 




where P and T are static values at the compressor face. Assuming P = P, 
combination of this equation with Eqn. (2) yields 

/ a T y-1 

^ ^ (P' ) Y = 1 + (e) Y 

= 1 + "'^Y e 2! 1 + 

The deviation of local speed of sound from average speed of sound is less 
than the deviation of local stagnation pressure from average stagnation 
pressure. Therefore a/a will be set equal to 1. 

Table II is a summary of these established relative magnitudes in 
terms of e. The velocities, Ug' and u^' , are placed at their absolute 
upper bounds and is set equal to for later verification. 



TABLE H RELATIVE ORDERS OF MAGNITUDE 
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Now Eqn. (7) can be written 
1 3P' 



yP' ar' 



( U '(0-,' - U '(0„ I = 

e z e j 



e e 



0 



.0 pO 
e e 



Vyip N au^' 



2tt cl 9t' 



(7a) 



Dropping terms of order e 
1 aP' 



yP' ar' 



- = 0 



w 



0 



— (-) 

yP' u^' \ ar' / 



( 10 ) 



Similarly for Eqn. (8) 
1 aP' 



yP'r' ae 



+ (u,' 



(A) 

2 r 



0 



0 



0 0 
e £ 






e e 



2n a at' 

0 2 
e e 



(8a) 



0)v 



yP 



M-) 

'u^'r'Vae / 



(11) 



The size of co^' was assumed to be of order e although at this point it 
is, in fact, unknown. In arriving at Eqns. (10) and (11), is 
multiplied by small terms {u^‘ and u^.'); therefore little if any error 
should be introduced here. Subsequent analysis establishes the size of 



“z' • 


From Eqn. 


(9) 




1 


aP' 






yP' 


az' "( 




- Uq'(o^ 


1 


1 


1 1 


1 1 


£0 


e 


£ £° 


£ £0 



Vjip N au^' 



2tt a at' 



I I 



0 

e e 



( 12 ) 
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Equation (12) contains most of the terms whose magnitudes were generally 
uncertain. All terms are of the order e so no consistent simplification 
can be made. This equation will not be used in the subsequent trans- 
formations from Pj to vorticity. 

C. APPLICATION OF RESULTS 

Equations (10) and (11) are the required relations for transforming 
total -pressure maps to circumferential and radial vorticity maps 
respectively. The procedure, as applied to the J85 map, is outlined 
below. 

The expressions 9P'/9r; in Eqn. (10) and 9P'/96 in Eqn. (11) are 
obtained in a similar fashion from Fig. A-2. For one particular radial 
line, values of P' are plotted against r', the slope at any point being 
the local value of 9P'/9r'. This is done for sufficient radials to 
adequately cover the entire map. Local values of P' can be taken direct- 
ly from the Pj map, and from these u^' can be calculated (assuming 
constant static pressure). Similarly in Eqn. (11), for one discrete 
radius, values of P' at a position, 0, are taken at intervals through 
360 degrees. These values are plotted, the local slopes being numerical- 
ly equal to 9P'/9^- Other values necessary to compute can be taken 
directly from the map or readily computed. See Appendix B for the 
calculations outlined above. 

The results of these calculations, depicted in Figs. 10 and 11 as 
o)p' and u)q' contour maps, present a two-dimensional picture of conditions 
at the compressor face. Since v*(vxw) =0 , it follows from Gauss' 
Divergence Theorem that 

// (u • r) dA = /// (v • u) dV = 0 
A V 
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Application of this theorem to the maps of and wq' contours yields 
a quantitative estimate of the axial vorticity, co^', which is reflected 
in Table II. These calculations are contained also in Appendix B, and 
the resulting map is given as Fig. 12. 

Radial and circumferential vorticity maps for th,e TF-30 total - 
pressure data in Fig. A-10 have also been generated. They are Figs. 

13 and 14. 
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FIG. 10 RADIAL VORTICITY CONTOURS, J85 
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FIG. II CIRCUMFERENTIAL VORTICITY CONTOURS, J85 
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FIG. 12 CONTOURS OF NET FLUX 

OF CU'z FOR J85 
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FIG. 13 RADIAL VORTICITY CONTOURS, TF-30 
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FIG. 14 CIRCUMFERENTIAL VORTICITY 

CONTOURS, TF - 30 
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IV. ANALYSIS OF VORTICITY MAPS 



A. RADIAL VORTICITY INLET EFFECTS 

In vector notation radial vorticity is aligned parallel to local 
blade spans of the compressor, + 03 ^ pointing outward radially according 
to the right-hand rule. The vorticity in Fig. 10 is concentrated in 
two areas, -to^ between 5 and 7 o'clock and at 8 o'clock. In 

analyzing the effects of these local concentrations it must be noted 
that four distinct situations can be represented by a vorticity map such 
as this. These are listed in Table III. Vortex-sheets are aligned par- 
allel with the flow and vortex filaments are convected with the flow. 



TABLE III 


Possible Combinations for Radial 


Inlet Vorticity 


Case 


Vortex Blades 


Flow 



I 


Sheet 


Stator 


Steady 


II 


Sheet 


Rotor 


Unsteady 


III 


Filament 


Stator 


Unsteady 


IV 


Filament 


Rotor 


Unsteady 



Case I 

A continuous sheet of vortex filaments passing through a stationary 
blade row (stator) is a steady-state problem. According to the Biot- 
Savart Law, velocities are induced in the fluid which cause changes in 
local blade incidences. Fig. 15. In a cascade r is the circulation 
about a blade. 1 
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( 



■■ = '(''ul - 



(13) 



and 



L = pV„r 



(14) 



Then 



aL Ap aV Ar 

— = — + — (15) 

L P r 




FIG. 15 RADIAL VORTEX -SHEET PASSING 
THROUGH A STATOR 



Therefore blades in region "A" where Ar/r is positive are more heavily 
loaded. 

The first stator row is behind a rotor so Figs. 10 and 13 cannot be 
applied here without knowing the influence of the rotor on tOj,'. This 
case could be applied directly to inlet guide-vanes however. 

Case II 

A steady vortex-sheet passing through a rotating blade row is an 
unsteady problem due to the relative motion between particle paths and 
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the rotor. Again, induced velocities cause variations in the absolute 
velocity, V-), which result in A3-] (relative inlet flow angle) and hence 
ar and aL. W is relative velocity in Fig. 16. As blade "A" passes 
through the vortex-sheet to position "B", it experiences a decrease in 
lift which is represented by a smaller ^lift* This results in shed 
circulation equal in direction to and a resultant decrease in 

Strength of the trailing vortex-sheet. Dynamic response of the blading 
is a factor which must be considered in this case. 

For a rotating cascade, circulation associated with blade lift is a 
function of the peripheral velocity, W^. For the circuit abed in Fig. 

16 this is 

'■lift = S(W„, - (16) 

The corresponding lift on this blade is 

L = prW^ = pS(W,i - W^^2)W„ (17) 




FIG. 16 RADIAL VORTEX -SHEET PASSING 
THROUGH A ROTOR 



48 




4 






% 



¥ 



For our sign convention, causes blade unloading while increases 
the local incidence. Comparison of circulation causing blade lift to 
that associated with inlet distortion (rdist^ should give 

a feel for the relative strength of this effect. For brevity, all mean 
quantities will be used and the dynamic response of the blading will be 
neglected. From data in Appendices C and D, for the J85 



= 650 ft/sec 



= 321 ft/sec 



Blade Gap (S^) = ft 



1 05 

Blade Chord (C) = ft 

12 




(16) 



^lift ■ ft^/sec (Relative) 



^dist 



(18) 




(D-1) 



From Fig. A-2 at 8 o'clock 



P 



avg 



From Fig. 10 at 8 o'clock 




+ .13 



Therefore 



and 



(Op = + 372 sec ^ (absolute) 

Tdist “ ft2/sec (absolute) 
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The ratio of circulations becomes 



lift 

dist 



30 

2.97 



10 



which is about 10 per cent. It must be remembered that this is occurring 
at blade passing frequency, and any induced stall would be dynamic stall. 



o, . r, . r /14,200 Rev\^31 Blades^ 

Blade Passing Frequency = 60~ Re\T7^ 



Blades 

ms. 



Case III 

Unsteady vortex filament motion through a stator is a more complex 
problem since time dependency of the fluid motion as well as the unsteady 
response of the blading must be considered. This situation can be 
pictured as the passage of a vortex filament or a local concentration of 
vorticity passing through a cascade. In a similar problem, Crimi and 
Reeves [Ref. 25] calculated sizable lift fluctuations as a helicopter 
rotor blade passed over a spanwise oriented vortex. Ward [Ref. 26] 
observed dynamic stall of a rotor blade which he asserted was caused by 
a radial vortex passing beneath the blade. Figure 17 depicts this 
situation in a cascade. 




FIG. 17 RADIAL VORTEX FILAMENTS 
PASSING THROUGH A STATOR 
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Case IV 



It follows that a vortex filament approaching a rotating blade row 
is also an unsteady problem. Figure 18 shows basically how W^J 1 is 
affected by Wp. 




FIG. 18 RADIAL VORTEX FILAMENT 
APPROACHING A ROTOR 

B. CIRCUMFERENTIAL VORTICITY INLET EFFECTS 

Circumferential vorticity taken vectorially is positive in the 
clockwise direction. Due to the signs in Eqn. 10, the positive gradients 
inFigs.ll and 14 represent -lo) 0 'l. The circumferential vorticity map 
for the J85 in Fig. 11 shows a strong concentration of vorticity at the 
outer radii extending through about 300 degrees. Figure 14 shows a 
similar characteristic for the TF-30 although it is weaker and less 
extensive. As in the case of tOp', four basic situations can exist for 
possible analysis. These are listed in Table IV. 
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TABLE IV 


Possible Combinations 


for Circumferential 


Inlet Vorticity 


Case 


Vortex 


Blades 


FI ow 


I 


Sheet 


Stator 


Steady 


II 


Sheet 


Rotor 


Steady 


III 


Filament 


Stator 


Unsteady 


IV 


Filament 


Rotor 


Unsteady 



Case I 

Circumferential vorticity is oriented perpendicular to the local 
blade spans of the compressor. A sheet of coq approaching a stationary 
blade row is pictured in Fig. 19. 




FIG. 19 CIRCUMFERENTIAL VORTEX- SHEET 
APPROACHING A STATOR 
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A 




I 



Induced velocities according to the Law of Biot-Savart cause Vg-j and 
hence a-j to change ahead of the cascade. This results in variable lift 
along the blade spans which produces shed vorticity downstream. Shed 
vorticity will be dealt with in more detail later when considering 
secondary flows. 

Case II 

A circumferential vortex sheet approaching a rotating blade row 
presents two problems for analysis. Relative motion between the particle 
paths and the rotor gives rise to a streamwise component of vorticity as 
well as a component perpendicular to the blade spans. Fig. 20. 




VORTEX- SHEET APPROACHING A ROTOR 

The concentration of co^' at the outer radius in Fig. 11 is essentially 
constant for about 300 degrees of angular displacement. The life time 
of one particular characteristic such as this can be as long as the time 
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for one revolution. The series of TF-30 Pj maps in Appendix A shows a 
life-time of close to 1/4 revolution. The significance of this is that 
one blade can experience the same velocity distortion for up to 360 
degrees of revolution which allows a steady-state analysis. Induced 
velocities caused by the normal component of this sheet are those shown 
in Fig. 19. Velocities induced by the streamwise component are pictured 
in Fig. 21. 




FIG. 21 STREAMWISE COMPONENT OF A 

CIRCUMFERENTIAL VORTEX- SHEET 
' AT A ROTOR 



The streamwise component causes the typical variations in lift along the 
blade span, but it also passes through the blade row and affects the 
exit flow angle. This will be discussed later. The normal component 
also causes a spanwise variation in lift as shown in Fig. 19. 

For ojq of small circumferential extent, the situation is unsteady 
due to relative motion between particle paths and the rotor. 
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I 



# 



Case III 



A circumferential vortex filament approaching a stator is an un- 
steady problem. The effect of the associated entry velocity perturbations 
shown in Fig. 22, a function of the dynamic response of the blading, would 
probably be negligible. However, as this filament passes through the 
blade passage, a streamwise component develops which can interact with 
existing secondary flows. This phenomenon will be discussed in the next 
section. 




FIG. 22 CIRCUMFERENTIAL VORTEX FILAMENT 
APPROACHING A STATOR 

Case IV 

A circumferential vortex filament approaching a rotor is different 
than the previous case only because a streamwise component of the filament 
is present ahead of the blades. Again unsteady response of the blades is 
a factor to consider. 
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c. PASSAGE AND EXIT EFFECTS 

Several authors have treated the problem of turning a non-uniform 
flow in a cascade of airfoils. Squire and Winter's method for calcu- 
lating secondary flows downstream of a cascade [Ref. 27] has been used 
extensively in the literature. Preston [Ref. 18], using a more simpli- 
fied approach, also treates the generation of streamwise vorticity which 
gives rise to secondary flow after a cascade. However, Preston assumes 
that in addition to the components of vorticity present in the original 
vortex filament, a change in circulation about the blade will cause 
circulation to be shed into the stream. Hawthorne [Ref. 16 and 17] 
attempts an analytical proof of the existence of these different com- 
ponents of vorticity in nonuniform flow, and he obtains reasonable 
agreement with experimental results. 

Consider an isolated vortex filament in a sheet normal to the flow 
direction as it passes through a cascade of stationary airfoils; see 
Fig. 23. Due to the difference in flow velocity over the suction and 
pressure surfaces of the blades, this filament will be turned and 
stretched as it passes through the cascade. Downstream of the cascade, 
segment ab' gives rise to a "distributed streamwise vorticity" through- 
out the passage exit plane. This description applies not only to a 
filament as shown here but also to a vortex-sheet. Squire and Winter s 
well known expression for this distributed vorticity is 



’dist 



■2e 



du 

dz 



(19) 



where i is streamwise vorticity, e is the deflection angle, and Z is 
the spanwise position. 
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FIG. 23 VORTEX FILAMENT PASSAGE 
THROUGH A CASCADE 



Segment b-b', part of the original filament a-b in Fig. 23, lies along 
the stagnation streamline and in considering the spanwise direction, this 
becomes a "vortex-sheet." This is called "trailing filament" vorticity. 

A third component of exit vorticity, "shed circulation" caused by vari- 
able spanwise lift due to the asymmetric inlet flow, also lies along the 
stagnation streamline and gives strength to the trailing vortex-sheet. 

This component is labeled c-d in Fig. 23. 

Preston's simple analytical expressions for these secondary components 
of circulation, as applied to Fig. 23, are as follows. 

Distributed Secondary Circulation : 



And replacing C 2 with Squire and Winter's expression for distributed 
secondary vorticity, Eqn. (19), yields 



dr-| = ^2 S cos (|.) dz 



dr 




(20) 
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Trailing Filament Circulation: 



‘"2 ^ 



(21 



Trailing Shed Circulation : 



'"'3 ^ ‘-r’ 



( 22 ) 



Trailing shed circulation is equal in magnitude and direction to dr 2 - 
The net circulation associated with the trailing vorticity in one blade 
passage and from one airfoil in the form of a "vortex sheet is the sum 
of these three components. 



From this exit secondary circulation the change in mean exit flow angle 
can be calculated. This total circulation is not uniformly distributed 
over the exit area [Ref. 18]. 

Investigations into changes in deflection due to cascade secondary 
flows have generally taken two lines of thought. Values of induced 
deflection at the lifting line v/ere derived from the induced velocities 
calculated solely from the strength of the vortex sheets by Carter 
[Ref. 28] and Hausmann [Ref. 29] to mention a few. These methods neglect 
the effects due to distributed passage vorticity and hence are unlikely 
to provide satisfactory agreement with experiment [Ref. 12]. 

Hawthorne [Ref. 30], Smith [Ref. 31], and others have considered 
secondary motions to be induced by the distributed vortices within the 
passages. It is assumed that the flow downstream behaves as if each 



dr = dr-| + 6T2 + dr3 



dr 




- 2 e S cos (|-) ] dz (23 
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blade passage were separated from its neighbor by the blade wakes, i.e. 
vortex sheets. The secondary velocities in the blade passages, v 
circumferentially and w spanwise, are defined by a stream function, ip, 
such that 

w = - ii V = -^ 

The stream function satisfied continuity for this steady secondary motion 
superposed on the main flow. 



9w ^ av 
9y 9z 



= 0 



Substituting the stream function into the vorticity equation (co = curl u) 
and setting it equal to Squire and Winter's expression for the distributed 
vorticity, Eqn. 19, yields Poisson's Equation. 



9y^ 9Z^ 




(24) 



Solutions to Eqn. (24) are tedious but have been accomplished by a series 
method [Refs. 27and 30], by relaxation techniques [Ref. 17], and by 
electrical analogs [Ref. 32]. The latter reference considers variable 
spanwise deflection which causes Eqn. (24) to become 



7^ip = 



= -2e 



du-] 

dz 



de 
^1 dl 



(25) 



An exit flow solution for a parabolic inlet profile similar to the case 
at hand is given by Honda [Ref. 33]. 
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Once secondary velocities, v and w, have been determined, local 
changes in the exit flow angle, Aa 2 > can be computed. V 2 is the main- 
stream exit velocity for uniform inlet flow with exit angle 02 * 




The exit flow angle for distorted inlet flow then becomes 

dist 

Comparison of ao with appropriate blade loading criteria [Ref. 10 
dist 

and 11] could explain the compressor's unstable operation. 

Another method of dealing with secondary velocities, v and w, ob- 
tained from solving Poisson's Equation was used by Hawthorne [Ref. 17]. 
Again the assumption of blade passages extended downstream by vortex- 
sheets trailing the blades is made. There are no normal velocities at 
the vortex-sheets. The strength of these sheets can be determined by 
calculating Aw, the change in spanwise secondary velocity across them. 

As mentioned before, these sheets of fluid leaving the trailing edges 
are vortex-sheets composed of trailing shed and trailing filament 
vorticity; see Fig. 24. If the trailing filament vorticity is estimated 
and subtracted from the total strength of the sheet, the trailing shed . 
circulation remains. From this can be calculated the variation in lift 
on the blade since shed circulation and blade lift are related by 

AL = Pur^ 

Hawthorne's expression for this trailing filament circulation is 

dr, = sin 2c, - sin 2„ 

\ dz / 1 2 cos ai 



- 2e 



COs(ag) 






(26) 
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Blade Trailing Edge 




FIG. 24 SECONDARY FLOW DOWNSTREAM OF A 
BLADE ROW [ adopted from Ref. I?] 



which is the same as Preston's formula given as Eqn. (21) except for 
different inlet/exit flow angles. Hawthorne's predicted shed circulation 
agreed quite well with an experimentally determined spanwise lift 
variation. 

D. CALCULATION OF SECONDARY CIRCULATION 

A calculation of these components of secondary circulation as caused 
by the circumferential vorticity in Fig. 11 was performed in order to 
quantitatively evaluate the effect of such distortion on a compressor. 

The analysis was made assuming steady-state conditions which means Case 
II of Table IV, a circumferential vortex-sheet passing through a rotor. 
The following terms require definition. 

V - Absolute Velocity 

W - Relative Velocity 

U - Rotor Peripheral Velocity I 
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0 ) 

C - 

fi - 



a 

3 - 






Subscripts! 
r,e,z- - 
s,n - 



0 - 

1 - 

2 - 



Absolute Fluid Vorticity 
Relative Fluid Vorticity 
Rotor Angular Speed 
Absolute Flow Angle 
Relative Flow Angle 
Velocity Vector 
Vorticity Vector 
see Fig. 25 

Absolute Coordinate Direction 
Natural or Streamline Coordinates 
s - Streamwise 
n - Normal to Streamlines 
Upstream of Compressor 
Compressor face (entry to blade row) 
Blade-row exit plane 




FIG. 25 NOMENCLATURE 



62 



Horlock's method [Ref. 15] of treating the passage of a vortex 
filament through a rotating blade-row results in the following vector 
transformation. Figure 26(a) is the velocity diagram ahead of the rotor 
which is convecting the vortex filament, wq (one of many in a sheet). 

The absolute velocity is axial so that Wg-j = Figure 26(c) is 

the inlet absolute vorticity combined with rotor rotation to give the 
inlet relative vorticity. Equations (27) and (28) are the corresponding 
analytical expressions for these components of relative vorticity in 
terms of Ug. 

= +2 n cos 3 -| - oj 0 sin 3 -| (27) 

= -2 n sin 3 -| - Wg cos 3 -| (28) 

The approach of a vortex filament to a rotating cascade and the resultant 
secondary relative vorticity at the passage exit are pictured in Fig. 27. 
This figure is similar to Fig. 23 except it applies to a rotor while Fig. 
23 is for a cascade. 

Combining inlet streamwise vorticity which passes through the cascade 
unaffected with the exit vorticity of Eqn. (19) yields the total dis- 
tributed vorticity in the blade passage exit plane. 

^s2 = ^sl ■ 2 e (29) 

Substituting Eqns. (27) and (28) 

?52 = (2fi cos 3-| - sin 3-| ) - 2e(-a)g cos 3-| - 2n sin 3-|) (30) 

Since this exit vorticity is in the direction of the flow, multiplying 
it by the blade passage exit area normal to the direction of flow yields 
exit distributed circulation. 
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(c) 



FIG. 26 VECTOR TRANSFORMATION FROM 

ABSOLUTE TO RELATIVE VORTICITY 
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FIG. 27 PASSAGE OF A VORTEX FILAMENT 
THROUGH A ROTATING CASCADE 

-= ^32 ^ 

When transformed from absolute to relative coordinates, Hawthorne's 
expression for trailing filament circulation given previously as Eqn. 
(26) becomes 



dr 



2 



dr 




sin 262 - sin 26-| 
cos e-| 



2e cos 3 




(32) 



Assuming two-dimensional flow, Hawthorne's equation for shed circulation 
is 



dr 



3 



dr 



-S cos 3-| ^.j(tan ^2 ~ ^l) 



(33) 



A simplified application of these equations to the circumferential 
vorticity and total-pressure maps for the J85 (Figs. 11 and A-2) was 
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made by assuming an average radial distribution of 3P'/9r' and P' to 
exist for a circumferential segment of the face. Equation D-2 was 
applied to radials at the three circumferential positions shown in 
Fig. 28. (Appendix D) 



U)c 



( -2854.9 \ 

\ py9.17P' - 8.17/ 



ap' 
ar ' 



sec 



-1 



(D-2) 



Now using Eqns. (27), (28), and (29), the various components of relative 
yorticity were calculated for a sufficient number of radial positions; 
see Appendix E. Figure 29 shows the relative sizes of these components 
for circumferential position I. 





FIG. 28 DATA POINTS FOR CALCULATION OF ^ 
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VORTICITY (sec"l) 




FIG. 29 VORTICITY COMPONENTS AT RADIAL I 
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Secondary circulations as a function of span at the blade row exit plane 
were calculated using Eqns. (31), (32), and (33); see Appendix E. The 
area under curve "A" in Fig. 30 represents distributed secondary circu- 
lation. The area under curve "B" is trailing filament circulation, and 
the area under "C" is shed circulation. These latter two are the bir- 
culations due to the vortex-sheet. The magnitudes of these components 
are found by integration over dr'. 

.9 

R } dr' = 32.1 ft^/sec 

.3 V dr / 

where R is the radius used for non-dimensionalizing. 

Concerning the results in Table V, distributed secondary circulation, 
r-| , represents a loss in kinetic energy. Blade-lift circulation is 
approximately equal to 30 ft^/sec; therefore the ratio of these two is 
about 1 to 1. In order to evaluate the actual effect r-| has on the 
compressor, a solution to Poisson's Equation, as described previously, 
would be necessary. From this solution the exact strength of the 
trailing vortex-sheet could be determined. Subtracting trailing fila- 
ment vorticity (or V 2 ) from the calculated strength of this vortex-sheet 
would yield the exact shed circulation. Equation (33) from which 
in Table V has been calculated was found by Hawthorne to yield values 
higher than those found experimentally. However, this analysis is 
essentially an order-of-magnitude check, and Eqn. (33) is adequate. 

Since shed circulation is related to the spanwise variation in blade 
lift by 

aL = p W„ 

Al _ ^shed 
L ' ^lift 



then 



(34) 
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I 





The circulation about one blade in the rotor under consideration for 
uniform inlet flow is approximately 

^lift = ft2/sec 

The average value of positive shed circulation over the three circumfer- 
ential data points considered, Radials I, II, and III, is 



5.0 + 3.9 + 15.1 
3 



8 ft^/sec 



So the variation in lift from Eqn. (34) is 



AL _ 8 

T ■ 30 



25 % 



which represents unloading at the center span of the blade. For the 
negative shed circulation at the blade tip 

-6.0 - 8.5 - 10.9 AL -8.8 



If the compressor blading is operating near stall it certainly seems 
plausible that a_ change of 25 to 30 per cent in lift could trigger stall. 
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Table V Magnitudes of Circulation Components 
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V. SUMMARY AND REVIEW 



Figures 31 and 32 are presented as a review and outline of the 
preceding theory for analyzing the effect of inlet flow distortion on 
an axial-flow turbomachine. Data from total -pressure maps, 8Pj/9r and 
9 Pt/ 9®» consitute the initial known quantities. All pertinent assump- 
tions and simplifying restrictions are included in the chart as they 
apply. The resulting determination of through a consistent order- 
of-raagnttude analysis implies complete definition of the velocity field 
ahead of the compressor with the exception of the ambiguity in sign of 
0 ) 2 . This is significant in that a three-dimensional flow field has been 
determined from a two-dimensional pressure map at the expense of intro- 
ducing an error of order Continuing in Fig. 32, the analysis next 
relates portions of several existing theories on secondary flov/s in 
cascades and turbomachines to the maps of radial and circumferential 
vortldty. This yields an estimate of the strength of the secondary 
flows caused by inlet flow distortion. The large value for aL/L which 
is calculated suggests that the linear cascade theories used are not 
adequate for the accuracy required in deriving a stall criterion. Some 
necessary considerations for further refinement of the theory which were 
not considered in this analysis are (1) blade twist, (2) three-dimensional 
effects and interaction between normal secondary flows with distortion 
induced flows, apd (3) changes in vorticity with passage through the 
blade row. 
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FLOW DIAGRAM FOR ANALYSIS; STEADY FLOW 




Fig. 32. Composite Theory for Analyzing Vorticity Maps 



VI. CONCLUSIONS 



This analysis is the first attempt at relating total -pressure 
distortion maps to the internal aerodynamics of a turbomachine. Separ 
rating the inlet distortion into components of radial and circumferential 
vorticity introduced flexibility in treating the problem as well as 
maintained a comprehensible physical picture of the aerodynamical 
situation. For the data analyzed herein, radial vorticity was found 
to be less significant to degradation of compressor performance than was 
circumferential vorticity since the radial vorticity was coupled to the 
dynamic response of the blading. This would require an unsteady 
analysis. The relative significance of and in this analysis does 
not necessarily apply in general. The effect of axial vorticity, 
must be included in further work. 

The size of the calculated shed circulation from a blade indicates 
distortion induced secondary flows are of sufficient size to seriously 
affect the compressor's operation. With further refinement to the 
present theory, secondary flows could be evaluated to a degree of 
accuracy capable of predicting stall through the mechanism of shed 
circulation or changes in exit flow angles. A solution to Poisson's 
Equation would provide the needed accuracy with which to use appropriate 
bladerloading criteria in defining inlet-distortion stall criteria. 
Although stall normally originates somewhere downstream of the first few 
stages, the distortion must be treated as it passes through each indi- 
vidual stage. The present theory deals only with the first blade row; 
however, it is one row closer to the origin of the stall and has elimi- 
nated a certain degree of uncertainty. This theory may be extended to 
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<101- 



I 



succeeding stages in the compressor by the technique illustrated in 
Fig. 26. 

Finally, in working toward a universal distortion factor, the 
effects of other types of inlet distortion, i.e. total -temperature and 
molecular weight, must be included. These extensions to the present 
theory should provide reliable and comprehensive stall -inducing 
distortion criteria. 
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APPENDIX A 



ORDER OF MAGNITUDE ANALYSIS 



Part I Preparation of Total -Pressure Data 

The following data and calculations apply to the experiments per- 
formed on a J85-GE-13 turbojet as described in Ref. 2, a = 0^ stall 
point. Some data were also obtained from Ref. 24. 

Subscripts: 0 - Free stream 

1 - Compressor face 

2 - Rotor exit, first stage 

T - Stagnation quantity 

Fjl = 1549 psf 

P^ = 1380 psf (Static) 

= 532° R 

M = 2.5 

00 

2 

Compressor Inlet flow area, A-| = 185 in 
Figure A-1 is a pressure contour map of Pji/^jq reproduced from 
Ref. 2. For our purposes, this must be converted to P' which is defined 
as Pji/Pji • 

Pressure Recovery of Inlet = .77 



^TO 

^T1 ^ Pfl ^TO 

^T1 ^TO ^T1 

^1 ^ '"ti 1 



(A-1) 
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FIG. A-1 PEAK INSTANTANEOUS PRESSURE 
CONTOURS, J85 (from Ref. 2 ) 
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Application of Eqn. (A-1) to Fig. A-1 yields the desired pressure 
contours. Fig. A-2. Figures A-3 through A-11 are instantaneous pressure 
contours (P‘) taken from experiments conducted on a TF-30 engine [Ref. 
23j. The radial lines at the centers represent compressor rotation. 

RPM = 10,000 

time between Figs. = / V ^ \ 

\ 10^ Rev A 32 Seg / 

= .1875 X 10'^ 

= . 2 ms . 

• Elapsed time (A-3 A-12) = 1.5 ms. 

Figure A-11 is the frame prior to engine surge. 



Fj = 1440 psf 

Estimated Inlet Values 

R = .4 






Fj = 1610 psf 
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Flow Into Paper 




FIG. A-2 TRANSFORMED PRESSURE 

CONTOURS 
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FIG. A-3 



FRAME ONE OF INSTANTANEOUS 
PRESSURE CONTOURS SEQUENCE, TF-30 
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FIG. A-4 



FRAME TWO OF INSTANTANEOUS 
PRESSURE CONTOURS, TF - 30 
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FIG. A- 5 



FRAME THREE OF INSTANTANEOUS 
PRESSURE CONTOURS SEQUENCE, TF-30 
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FIG. A-6 



FRAME FOUR OF INSTANTANEOUS 
PRESSURE CONTOURS SEQUENCE, TF-30 
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FIG. A-7 FRAME FIVE OF INSTANTANEOUS 

PRESSURE CONTOURS SEQUENCE, TF - 30 
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FIG. A-8 



FRAME SIX OF INSTANTANEOUS 
PRESSURE CONTOURS SEQUENCE, TF-30 
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FIG. A-9 FRAME SEVEN OF INSTANTANEOUS 

PRESSURE CONTOURS SEQUENCE . TF- 30 
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FIG. A-IO 



FRAME EIGHT OF INSTANTANEOUS 
PRESSURE CONTOURS SEQUENCE, TF-30 



88 






FIG. A-ll 



FRAME NINE OF INSTANTANEOUS 
PRESSURE CONTOURS SEQUENCE, TF -30 
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Part II Pressure Fluctuations in the Z Direction 



In the case of steady flow , aPj is present only when dissipation 
occurs. Some causes of such dissipation in an inlet are shock waves, 
flow separation, obstructions, and boundary layers. These examples 
are illustrated in Fig. A-12. 



Boundary 

Layer 




Flow 

Seporotion 






(b) Wall Effects 









oundary Layer 



Strut 




FIG. A-12 SOURCES OF DISSIPATION 

WITHIN INLETS 
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The stagnation pressure gradient is present only in the vicinity of dis- 
sipation. Therefore in steady flow, aP'/Sz' is zero at the compressor 
face if the dissipation is located upstream as shown in Fig. A-13. 




FIG. A-13 PRESSURE HISTORY FOR 
STEADY INLET FLOW 

For unsteady flow , pressure gradients are again caused by dissipation. 
If there were a standing shock wave in the inlet, unsteady flow could 
cause this to oscillate in the axial direction. This is called "shock 
wave twitter" and can be self excited in oscillation. This is depicted 
in Fig. A-14. 




FIG. A-14 SHOCK- WAVE TWITTER 



91 



The changing geometry seen by this oscillating shock-wave causes pres- 
sure pulses to be generated as shown in Fig. A-15, A/A* is the area 
ratio as a function of x. The second quadrant is local Mach number as 
a function of area ratio, A/A*. The third quadrant is the stagnation 
pressure ratio for a shock at M. Finally, the fourth quadrant yields 
the derivative of Pj with respect to x. 




FIG. A-15 TOTAL PRESSURE FLUCTUATIONS 
FROM OSCILLATING SHOCK- WAVES 



The resulting aP^/ax, a function of diffuser design, can be evaluated 
at the compressor face by assuming an acoustic analysis with no wave 
reflections at the compressor. Figure A-16 is a plot of the pressure 
pulses (waves) in the x-t plane. The slope of the characteristics is 
equal to the reciprocal wave speed; curvature is due to the variable 
area at the center-body. 
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FIG. A-16 PRESSURE WAVES IN THE x-f PLANE 
AS A RESULT OF SHOCK-WAVE TWITTER 



Consider the three points labeled 1 , 2, and 3 in Fig. A-16. 



aP' ? 2 ' - P^' 
3Z' - 22* - ' 

Since P-j' = P 2 ' 



?P' P 3 ' - P^ ■ ? 2 ' - Pi ' 

^ ‘3' ■ ‘1' » 3 ' - t,' 



equating 
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22* -Z^' (Z2 - Z^)(2/D) 1 (z2-^i) 

ts'-ti' (t3 - ti)(2?/D) I 



1 

_ (a + u) = a' + u' 
a 



therefore 



aP' 1 

az' (a' + u') 




and 




1 

(a + u) 




(A-2) 



The usefulness of Eqn. (A-2) can be demonstrated by application to a 
simple inlet/engine configuration with an oscillating shock wave present; 
see Fig. A-17. For brevity the following assumptions are made 
(1 . ) a = Constant 

(2.) Wave impedance of compressor matches wave impedance 
of duct 

(3.) U = constant 
(4.) No wave dispersion 



L H 




FIG. A-17 NOMENCLATURE 
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The shock wave, oscillating within tN, produces pressure waves which 
are represented as a sine wave function. 



X - Xq = N sin (tot) 



8t 



to N cos (tot) 



within an 



aPj 9X 8P 

• 1 

at " at ax 

3 P 

= to N cos (tot) ( I) 

\ax L 
^0 

Using a wave transit time equal to L/(a+u) with a = wave speed relative 
to the gas, the unsteady pressure at the compressor face becomes 



at 



N cos to(t - 



Xq+L 






Combining this with Eqn. (A-2) yields 



3 Pi 



3z 



to 



N COS to(t - -i:^) / aP- 
a+u / T 



face _ 



u + a 



ax 



Non-dimensionalized the preceding equation becomes 

to / 1 ^ ^ ^ 

3P' (to/t2)(t2N/a) cos [ ^ a* + u^ ^ V 

az' Li' 

where Q is rotor angular speed. 

This is one method of relating pressure fluctuations at the compressor 
face (aP'/3z‘) with the pressure waves generated by inlet/shock-wave 
interaction. 
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APPENDIX B 



CONVERSION OF Pj MAP TO VORTICITY - J85 



Pressure data taken directly from Fig. B-1 are listed in Table B-I. 
From these data the curves in Fig. B-2 are generated. Slopes taken 
from these curves, which are tabulated in Table B-II, provide a means 
for constructing the circumferential vorticity map in Fig. 11. A 
similar procedure is performed in generating the radial vorticity map 
in Fig. 10. Data are contained in Tables B-III and B-IV and Fig. B-3. 
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FIG. B-l PRESSURE DATA AND SCALES FOR J85 
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Table B-I Pressure Data Taken in Radial Direction from J85 



Total-Pressure Map 



0' 

r' 


P' 


30 

r' 


0 

P! 


60' 

r' 


P' 


90 

r' 


O 

p. 


120® 

r' P' 


150° 

r' P' 


.4 


.91 


.47 


.91 


.46 


.91 


.44 


.91 


.71 


.91 


.47 


.84 


.68 


.97 


.7 


.97 


.75 


.97 


.62 


.91 


.75 


.97 


.55 


.84 


.72 


1.04 


.74 


1.04 


.79 


1.04 


.72 


.91 


.79 


1.04 


.68 


.91 


.78 


1.10 


.79 


1.10 


.82 


1.10 


.75 


.97 


.82 


1.10 


.75 


.97 


.81 


1.17 


.81 


1.17 


.84 


1.17 


.79 


1.04 


.84 


1.17 


.78 


1.04 


.92 


1.17 


.92 


1.17 


.92 


1.17 


.82 


1.10 


.92 


1.17 


.80 


1.10 


.95 


1.04 


.95 


1.04 


.95 


1.04 


.85 


1.17 


.95 


1.04 


.82 


1.17 














.92 


1.17 






.92 


1.17 














.95 


1.04 






.95 


1.04 
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FIG. B-2 
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FIG. B-2 (cont.) 
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FIG. B-2 (cont) 
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Table B-II Pressure Gradients Obtained from Slopes in Fig. B-2 



0° 

r' 3P72r' 


30 

r' 


0 

9P73r' 


60 

r' 


O 

3P79r' 


90° 

r' 3P79r' 


120° 

r' 3P79r' 


150° 

r' 3P79r' 


.4 


0 


.4 


.13 


.4 


.07 


.32 


.045 


.4 


.03 


.3 


-.22 


.5 


.09 


.5 


.16 


.5 


.12 


.38 


.17 


.5 


.06 


.3 


-.22 


.6 


.28 


.6 


.26 


.6 


.135 


.44 


.26 


.6 


.1 


.5 


-.22 


.7 


.9 


.7 


.5 


.7 


.25 


.52 


0 


.68 


.33 


.51 


0 


.76 


1.8 


.76 


1.4 


.74 


.56 


.6 


-.26 


.72 


1.08 


.54 


.3 


00 


2.6 


.79 


5 


.8 


2.3 


.62 


-.27 


.8 


3 


.6 


.53 


.80 


1.6 


.8 


6 


.82 


4.5 


. 66 


0 


.82 


2.5 


.7 


.73 


.82 


.83 


.82 


1.1 


.86 


.55 


.7 


.29 


.84 


1.08 


.76 


1.5 


00 


.37 


.86 


0 


-.89 


0 


.72 


1.2 


.86 


.4 


.79 


4.0 


00 


0 


.9 


-.72 


.92 


-.6 


.75 


2.0 


.89 


0 


.8 


7 


.9 


-.6 


.94 


-2.7 


.94 


-2.7 


.8 


2.0 


.9 


-.26 


8.2 


1.3 


.94 


-2.7 


.96 


-3 


.96 


-3.0 


.83 


6.0 


.92 


-.6 


.86 


0 


.96 


-3.0 










.85 


.9 


.94 


-2.7 


.9 


-.72 














.87 


.2 


.96 


-3.0 


.94 


-2.7 
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Table B-II (Continued) 
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Table B-III Pressure Data Taken in Circumferential Direction from J85 

Total-Pressure Map 



r' = . 


38 


• r' = 


.46 


r' = . 


54 


r' = 


.61 


r' = 


.69 


r' = 


.77 


0 

(rad) 


P' 


0 

(rad) 


P' 


0 

(rad) 


P' 


0 

(rad) 


P' 


0 

(rad) 


P' 


0 

(rad) 


P' 


4.26 


.91 


4.19 


.91 


1.81 


.91 


1.69 


.91 


0 


.97 


0 


1.04 


6.17 


.91 


4.48 


.97 


2.35 


.84 


2.41 


.84 


1.25 


.91 


.872 


.97 






5.23 


.97 


3.14 


.84 


3.19 


.84 


2.79 


.91 


2.09 


.97 






.174 


.91 


4.10 


.91 


3.98 


.91 


3.84 


.91 


2.96 


.97 










4.41 


.97 


4.36 


.97 


4.31 


.97 


3.40 


.91 










5.48 


.97 


5.72 


.97 


4.59 


1.04 


3.87 


.91 


















5.13 


1.04 


4.27 


.97 






















4.48 


1.04 





r' = 


.85 




r' = 


.93 


0 

(rad) 


P' 


0 

(rad) 


P' 


0 

(rad) 


P' 


.174 


1.17 


3.70 


v91 


3.03 


1.17 


.785 


1.10 


3.92 


.91 


3.19 


1.10 


1.39 


1.04 


4.22 


.97 


3.37 


1.04 


2.0 


1.10 


4.39 


1.04 


3.58 


.97 


2.26 


1.17 


4.52 


1.10 


4.10 


.97 


2.61 


1.17 


4.71 


1.17 


4.27 


1.04 


2.87 


1.10 


4.88 


1.17 


4.45 


1.10 


3.14 


1.04 


5.23 


1.17 


4.57 


1.17 


3.40 


.97 


6.11 


1.17 
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FIG. B-3 (conf.) 



Table B-IV Pressure Gradients Obtained from Slopes in Fig. B-3 
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Application of Gauss' Divergence Theorem is performed in the follow- 
ing manner. A grid such as is shown in Fig. B-4 is placed over a 
vorticity map. 




FIG. B-4 METHOD OF DIVERGENCE 

APPLICATION 



The net passing through region abed is equal to 

[(-.13)(^) + (fcJl)(^)] - [(J|^)(|cd) + (S^Xfed)] = 

Net flux of Wg' is calculated in a similar manner using the same grid 
over an Wg' map. The sum of flux of and Wg' is then equal to 
“^'(ab-cd) so that 

to,' (ab’cd) 

(O ' = 

Area (abed) 

This procedure is performed on each of the 36 segments so that a complete 
toz' map can be drawn. These data are contained in Table B-V. If the sum 
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of net flux of wp' and Wg' is negative, then the net flux of must 
be positive. There is an ambiguity of algebraic sign which cannot be 
resolved without additional data or information. One evaluates . 
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Table B-V Net Fluxes of Vorticity in Three Coordinate Directions 
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APPENDIX C 



ESTIMATION OF SOME J85 PARAMETERS 



Subscripts: T - Stagnation Quantity 

0 - Freestream ahead of Compressor 

1 - Compressor Face 

2 - Rotor Exit, first stage 
m - Mean radius 

Data obtained from Ref, 24: 

Pjl = 1549 psf 
= 1380 psf 
Tji = 532‘’R 

Data obtained from Ref. 2 and Ref. 34: 



« - ’’t 

® " 2116 psf 

N* = 16,500 RPM 



(Corrected Total -Temperature) 



(Corrected Total-Pressure) 



(Rated Engine Speed) 



Wcorr 



6 



15 KG/sec = 33.5 Ibs/sec (Corrected Mass Flow Rate) 



Percentage Corrected Engine Speed = x 100 = 86.9 



D-| = 16.1 inches 

^1 hub = inches 

A^ = 155 in2 



N = 



(Compressor Face Outside Diameter) 



31 
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(Entrance Flow Area) 

(Number of blades in first row) 






■ .F 



♦ t 






F[G. C“1 BLADE PASSAGE DIMENSIONS 



Calculated Parameters 




(1 + III m2) 
Y 




1549 

1380 



1.125 



(1 +,2M2) 



= .414 




’Ti 



1 + ^ 




514°R 







il 



I 



a 



49.02 



= yRTi 

F = 

¥ = 1110 ft/sec 



= M a 

= {.414)(1110 ft/sec) 

= 460 ft/sec 

F, = ^ 

rt 

1 - = 1380 Ibs/ft^ 

(1715.6 ) (514°R) 

^ = .001565 s1uqs/ft^ 



As a check on these figures, another method may be used. 



W 



1549 

l^cor g ^ (33.5 Ibs/sec) ^ ?il6 ^ 

yr v^532/518 



W = 24.15 Ibs/sec 



W = pA^U^ 




(24.15 Ibs/sec) (144 in^/ft^) 

(.001565 llMl ) (155 in^) (32.2 ) 

ft^ slug 

= 462 ft/sec 
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Estimation of 




FIG. C-2 VELOCITIES AHEAD 

OF ROTOR 



W Relative Velocity 
V -> Absolute Velocity 



6 , = 



tan 



-1 



(- 



D • RPM 



Dj^ - .875 ft 

Ang. Speed = j 86.9)(1650 ^ _ 14^200 RPM 

Vfj/518.7 X 100 



W 



ul 



Ti (.875 Jl- ) ( 14200 revN 
rev 60 sec^ 



W^ji = 650 ft/sec Mean Radius 



^1 



tan"1 



650 > 
460 ^ 



3-| = 55® Mean Radius 
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Estimation of ^2 

Individual stage characteristics and loss data for the J85 are not 

available. The following data are from Ref. 34: 

Compressor of Free-Vortex Design 

Lieblein Diffusion Parameter, D < .5 

Solidity (a) = 1.14 'j 

V First Stage Rotor 
Axial Velocity Ratio = .96 J 



D = 0 - ^) 

Wl 2YW^ 

W 2 = Va2/cos ^2 
= 442/ cos ^2 

Wu 2 = Vg 2 tan 62 
= 442 tan ^2 

' W-] = Vg/cos e-] 



(Ref. 10) 

.96 V^-j/cos ^2 
546 

.96 Vgi tan ^>2 
802 ft/sec 



442/cos &2 ) + ■ 442 tan 32 ) 

802 (2)(1.14)(802) 



32 = 36® minimum 
Wu 2 = 320 ft/sec 

Since the loading limit is D < .5, 36° can be considered a minimum for 
32 » and the maximum turning angle becomes 
43 = 3^-32 = 55° - 36° 

= 19° 

e = .332 radians 
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APPENDIX D 



CALCULATION OF ABSOLUTE VORTICITY 



Radial Vorticity: 



b) 



1 



yP'Uz' 



V r' 39 / 



(11) 



Y = 1.4 



P' = 



r\_ 

'ti 



Uz' 



u 



U: 



z _ 

a Uz/R 



= M ^ 
U, 



. ff /(PrPsmV' = M flLS. 



1380 



\l/2 



\(Pt-Ps)(2/p)/ V549 - 1380 



= r(' 



Jn - 8.17 

169 



) 



1/2 






= Mf "T ^T 



Pj 169 



\l/2 _/ P. \ 1/2 

8.17j = M ^9.17 pi - 8.17j 



U ' = My9/l7P' - 8.17 



03| 



I - 



2a 



Substituting into Eqn. (11) 



j = 2a / 1^ 3P' \ ft/sec 

^ 1.4 D P' M V"9.17P'-8.17 ( r' 39 j ft 
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\ 

( 






Recal 1 



a 

D 

M 



1110 ft/sec 
16.1 



12 

.414 



ft 



0 ) 



2854.9 



P' y9.17P' - 8.16 



n 5P' \ 
V r' 96 / 



sec 



-1 



(D-1) 



Circumferential Vorticity 



U)0 



-2854.9 / 9P' \ 

P' v'9.17P' - 8.16 \ ar' j 



sec 



(D-2) 



Equation (D-2) is used to convert the contours in Fig. 11 to 
absolute vorticity for subsequent calculation of secondary circulation. 
Tabulated results are contained in Table D-I. 
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Table D-I Absolute Vorticity for Three Circumferential Positions 



on cog' Map 







Radial 


I 






Radial 


II 




r' / 


r (in) 


P' 


3p'/ar' 


o)g (sec'^) 


P 


1 


3P'/3r' 


0 


(sec"l ) 


.3 


2.43 


.91 


.35 


-2554.9 


.9 


.05 


-520 


.325 


2.63 


.93 






-1770.8 






.05 






.35 


2.83 


.94 






-1567.6 






.05 






.375 


3.03 


.95 






-1416.3 






.05 






.4 


3.23 


.96 






-1297.8 






.05 






.425 


3.43 


.97 






-1201.6 






.06 


-624.1 


.45 


3.63 


.98 






-1121.4 


.! 


9 


.075 


-780.1 


.475 


3.83 


.99 






-1053.2 


.91 


.09 


-656.9 


.5 


4.03 


.99 






-1053.2 


.91 


.1 


-729.9 


.525 


4.23 


1.0 






-994.2 


.91 


.11 


-802.9 


.55 


4.43 


1.0 






-994.2 


.92 


.13 


-767.3 


.575 


4.63 


1.01 






-942.5 


.92 


.14 


-826.3 


.6 


4.83 


1.01 






-942.5 


.93 


.15 


-758.9 


.625 


5.03 


1.02 






-896.7 


.93 


.17 


-860.1 


.65 


5.23 


1.03 






-855.7 


.94 


.20 


-895.7 


.675 


5.43 


1.04 


.35 


-818.8 


.94 


.23 


-1030.1 


.7 


5.63 


1.05 


.75 


-1682.7 


.95 


.25 


-1011.6 


.725 


5.83 


1.06 


1.0 


-2156.2 


.97 


.5 


-1716.6 


.75 


6.04 


1.08 


1.: 




-2402.2 


1.0 


1.0 


-2840.7 


.775 


6.24 


1.10 


1.: 


? 


-2243.5 


1.04 


1.5 


-3509.2 


.8 


6.44 


1.17 


1.0 


-1522.4 


1.10 


2.0 


-3739.2 


.825 


6.64 


1.18 


.75 


-1112.4 


1.17 


1.5 


-2283.6 


.85 


6.84 


1.19 


0 




0 


1.18 


.5 


-741.6 


.875 


7.04 


1.20 


-.25 


352.6 


1.20 


0 




0 


.9 


7.24 


T.18 


-.75 


1112.4 


1.18 


-.5 


+741.6 


.925 


7.44 


1.17 


-1.0 


1522.4 


1.17 


-1.0 


1522.4 


.95 


7.65 


1.14 


-1.! 


5 


2480.2 


1.1 


-3.0 


5608.9 


.975 


7.85 


1.1 


-3.0 


5608.9 


1.05 


-3.2 


7179.8 


1.0 


8.05 





















118 



Table D-I (continued) 



r' / 


r (in) 


P' 


.3 


2.43 




1 


.325 


2.63 






.35 


2.83 






.375 


3.03 






.4 


3.23 






.425 


3.43 






.45 


3.63 






.475 


3.83 






.5 


4.03 






.525 


4.23 






.55 


4.43 






.575 


4.63 






.6 


4.83 






.625 


5.03 






.65 


5.23 


.90 


.675 


5.43 


.90 


.7 


5.63 


.91 


.725 


5.83 


.94 


.75 


6.04 


.97 


.775 


6.24 


1.0 


.8 


6.44 


1.04 


.825 


6.64 


1.17 


.85 


6.84 


1.18 


.875 


7.04 


1.2 


.9 


7.24 


1.18 


.925 


7.44 


1.17 


.95 


7.65 


1.04 


.975 


7.85 


1.0 


1.0 


8.05 







Radial III 

3P'/9r' o)„ (sec"^) 

0 



-.2 +2080.3 



-.2 


+2080.3 


0 


0 


+.3 


-3120.5 


.4 


-4160.7 


.5 


-5200.8 


.5 


-5200.8 


.5 




.5 




.5 


-3649.9 


1.0 


-4478.9 


1.5 


-5149.8 


2.0 


-5681.4 


2.3 


-5380.8 


1.7 


-2588.0 


0 


0 


-.3 


+423.2 


-.75 


+1112.4 


-1.5 


+2283.6 


-3.0 


+7018.4 


-3.2 


+9090.3 



APPENDIX E 



CALCULATION OF CIRCULATION FROM wq MAP FOR 085 



;51 = + 2fi cos 3^ - 0 )^ sin 3-| 

= + 2 (236.6 cos 55° - cog sin 55° 

^51 = + 271.4 - .8191 cog 



= - 2fi sin 3-| + u)g cos 3-| 



= 2 (236.6 sin 55° + (cos 55°) 



Sec 



= - 387.6 - .573 co. 



?32 = + ^ 



SI - 2e^nl 



= (271.4 - .8191 ug) + 2(.3316 rad)(387.6 + .573 cog) 



?S2 = 528.4 - .439 



dr 



= ?52 S cos ^2 



S = 



_ 1.095 in 



12 in/ft 



= .09125 ft 



= 36' 



dr- 



dr 



,0738 ^52 ft/sec 



(28) 

(E-1) 

(29) 

(E-2) 

(29) 

(E-3) 

(31) 



(E-4) 



120 



(32) 



dr2 

W~ 




sin 2$2 - sin 23-| 
cos 



2e cos 3 



2) 



= 55 ® 
e = .3316 



dF2 
dr ' 



.0480 ft/sec 



a 



(E-5) 




S cos g-j (tan B 2 " tan B-j) 



(33) 




.0336 



(E-6) 



Tabulated results for Eqns. (E-1) through E-6) using previously 
determined data are contained in Tables E-I, E-II, and E-III. 
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Table E-I Relative Vorticity and Circulation Components at Radial I 



r' 


0)0 

{sec~' ) 


?S1 


^nl 


^S2t 


dr.j/dr 


dF2/dr 


dr^/dr 




(sec" ' ) 


(sec"l ) 


(sec"l ) 


(ft/sec) 


(ft/sec) 


(ft/sec) 


.3 


-2554.9 


-H2364.1 


+1076.3 


+1650 


+121.7 


51.2 


+36 . 3 


.325 


-1770.8 


+1721.8 


+627.0 


+1305.7 


+96.3 


29.8 


+21.1 


.35 


-1567.6 


+1555.4 


+510.6 


+1216.5 


+89.7 


24.5 


+17.2 


.375 


-1416.3 


+1431.4 


+423.9 


+1150.1 


+84.8 


20.1 


+14.3 


.4 


-1297.8 


+1334.4 


+356 . 0 


+1098.1 


+81.0 


16.9 


+12.0 




-1201.6 


+1255.6 


+300.9 


+1055.9 


+77.9 


14.3 


+10.1 


.45 


-1121.4 


+1189.9 


+254.9 


+1020.6 


+75.3 


12.1 


+8.5 




-1053.2 


+1134.0 


+215.8 


+990.7 


+73.1 


10.3 


+7.3 


.5 


-1053.2 


+1134.0 


+215.8 


+990.7 


+73.1 


10.3 


+7.3 




-994.2 


+1085.7 


+182.0 


+964.8 


+71.2 


8.6 


+6.14 


.55 


-994.2 


+1085,7 


+182.0 


+964.8 


+71.2 


8.6 


+6.14 




-942.5 


+1043.4 


+152.4 


+942.1 


+69.5 


7.2 


+5.12 


.6 


-942.5 


+1043.4 


+152.4 


+942.1 


+69.5 


7.2 


+5.12 




-896.7 


+1005.8 


+126.2 


+922.0 


+68.0 


6.0 


+4.2 


.65 


-855.7 


+972.3 


+102.7 


+904.0 


+66.7 


4.8 


+3.5 




-818.8 


+942.0 


+81.5 


+887.8 


+65 . 5 


3.9 


+2.75 


.7 


-1682.7 


+1649.6 


.+576.5 


+1267.1 


+93.5 


27.4 


+19.4 




-2156.2 


+2037.5 


+847.9 


+1474.9 


+108.8 


40.3 


+28.7 


.75 


-2402.2 


+2239.0 


+988.8 


+1582.9 


+116.8 


47.0 


+33.3 




-2243.5 


+2109.0 


+897.9 


+1513.2 


+111.6 


42.6 


+30.3 


.8 


-1522.4 


+1518.3 


+484.7 


+1196.7 


+88.3 


23.0 


+16.3 




-1112.4 


+1182.5 


+249.8 


+1016.7 


+75.0 


11.8 


+8.4 


.85 


0 


+271.4 


-387.6 


+528.4 


+38.9 


-18.4 


-13.0 




352.6 


-17.4 


-589.6 


+373.6 


+27.5 


-28.0 


-19.8 


.9 


1112.4 


-639.7 


-1025.0 


+40.0 


+2.9 


-48.8 


-34.5 




1522.4 


-975.5 


-1259.9 


-139.9 


-10.3 


-60.0 


-42.5 


.-95 


2480.2 


-1760.1 


-1808.7 


-560.4 


-41.3 


-86.0 


-61 .0 


1.0 


5608.9 


-4322.8 


-3601.4 


-1933.9 


-142.7 


-171.5 


-121.0 
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Table E-II Relative Vorticity and Circulation Components at Radial II 



r' 


(sec 


4S1 

(sec"^ ) 


^nl 

(sec'l ) 


^S2 

j[sec~l ) 


dr-j/dr 

(ft/sec) 


.3 


' -520 


+697.3 


-89.6 


+756.6 


+55.8 




-520 


















.35 


-520 




















-520 


















.4 


-520 




















-624.1 


+782.3 


-29.9 


+802.3 


+59.2 


.45 


-780.1 


+910.3 


+59.3 


+870.8 


+64.2 




-656.9 


+809.4 


-11.1 


+816.7 


+60.2 


.5 


-729.2 


+869.2 


+30.6 


+848.8 


+62.6 




-802.9 


+929.0 


+72.4 


+880.8 


+65.0 


.55 


-767.3 


+899.8 


+52.0 


+865.2 


+63.8 




-826.3 


+948.2 


+85.8 


+891.1 


+65.7 


.6 


-758.9 


+893.0 


+47.2 


+861.5 


+63.5 




-860.1 


+975.9 


+105.2 


+905.9 


+66.8 


.65 


-895.7 


+1005.0 


+125.6 


+921.5 


+68.0 




-1030.1 


+1115.1 


+202.6 


+980.6 


+72.3 


.7 


-1011.6 


+1100.0 


+192.0 


+972.4 


+71.7 




-1716.1 


+1677.0 


+595.7 


+1281.7 


+94.5 


.75 


-2840.7 


+2598.2 


+1240.1 


+1775.4 


+131.0 




-3509.2 


+3145.7 


+1623.1 


+2068.9 


+152.6 


.8 


-3738.6 


+3334.5 


+1755.1 


+2170.0 


+160.1 




-2283.6 


+2141.8 


+920.9 


+1530.9 


+112.9 


.85 


-741.6 


+878.8 


+37.3 


+853.9 


+63.0 




0 


+271.4 


-387.6 


+528.4 


+38.9 


.9 


+741.6 


-336.0 


-812.5 


+202.8 


+14.9 




+1522.4 


-975.5 


-1259.9 


-139.9 


-10.3 


.95 


+5608.9 


-4322.8 


-3601.4 


-1933.9 


-142.7 




+7179.8 


-5609.5 


-4501 .6 


-2623.5 


-193.6 


1.0 





















dF2/dr 

(ft/sec) 

-4.26 



-1.4 

+ 2.8 

-.5 

+1.4 

+3.4 

+2.4 

+4.0 

+ 2.2 

+5.0 

+5.9 

+9.6 

+9.1 

+28.3 

+59.0 

+77.2 

+83.5 

+43.8 

+1.7 

-18.4 

-38.6 

-59.9 

-171.4 

-214.2 
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dr^/dr 

(ft/sec) 



-3.0 



- 1 .0 
+1 .9 
-.3 
+ 1.0 
+2.4 
+1.7 
+ 2.8 
+ 1.6 
+3.5 
+4.2 
+ 6.8 
+6.4 
+ 20.0 
+41.7 
+54.6 
+59.1 
+31.0 
+ 1.2 
-13.0 
-27.3 
-42.4 
-121.3 
-151.7 



Table E-III Relative Vorticity and Circulation Components at Radial III 



r' 



.3 



.35 

.4 



“ft 

(sec“^ ) 




^nl, 

(sec"^ ) 



^S2 dr-|/dr 
(sec"l) (ft/sec) 



dr2/dr 

(ft/sec) 



+ 2080.3 



- 1432.5 



- 1579.6 



- 384.8 



- 28,4 



- 75.1 



.45 



.5 

.55 

.6 

.65 

.7 

.75 



.8 



.85 



.9 

.95 



1.0 



0 

- 3120.5 

- 4160.7 

- 5200.8 



- 3649.9 

- 4478.9 

- 5149.8 

- 5681.4 

- 5380.8 

- 2588.0 

0 

+ 423.2 

+ 1112.4 

+ 2283.6 

+ 7018.4 

+ 9090.3 



+ 271.4 

+ 2827.4 

+ 3679.4 

+ 4531.3 



+ 3261.0 
+ 3940.0 
+ 4489 . 6 
+ 4925.0 
+ 4678.8 
+ 2391.2 
+ 271.4 
- 75.2 
- 639.7 
- 1599.0 
- 5477.3 
- 7174.4 



- 387.6 

+ 1400.4 

+ 1996.4 

+ 2592.4 



+ 1703.7 

+ 2178.8 

+ 2563.2 

+ 2867.8 

+ 2695.5 

+ 1095.3 

- 387.6 

- 630.0 

- 1025.0 

- 1696.1 

- 4408.1 

- 5596.3 



+ 528.4 

+ 1898.2 

+ 2354.9 

+ 2811.5 



+ 2130.7 

+ 2494.6 

+ 2789.1 

+ 3022.5 

+ 2800.5 

+ 1664.5 

+ 528.4 

+ 342.6 

+ 40.0 

- 474.1 

- 2552.6 

- 3462.2 



+ 38.9 

+ 140.0 

+ 173.7 

+ 207.4 



+ 157.2 

+ 184.1 

+ 205.8 

+ 223.0 

+ 213.3 

+ 122.8 

+ 38.9 

+ 25.2 

+ 2.9 

- 34.9 

- 188.3 

- 255.5 



- 18.4 

+ 66.6 

+ 95.0 

+ 123.3 



+ 81.0 

+ 103.7 

+ 122.0 

+ 136.5 

+ 128.3 

+ 52.1 

- 18.4 

- 29.9 

- 48.7 

- 80.7 

- 209.8 

- 266.3 



dF3/dr 

(ft/sec) 



- 53.2 



- 13.0 

+ 47.1 

+ 67.2 

+ 87.3 



+ 57.4 
+ 73.4 
+ 86.3 
+ 96.6 
+ 90.8 
+ 36 . 9 
- 13.0 
- 21.3 
- 34.5 
- 57.1 
- 148.5 
- 188.5 
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